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We discuss mechanisms of self-organization processes in two-level solid-state class-B laser system. 
The model is considered under assumptions of influence of nonlinear absorber and external force, 
separately. It was found that self-organization occurs through the Hopf bifurcation and results to 
a stable pulse radiation. Analysis is performed according to the Floquet exponent investigation. 
It was found that influence of the nonlinear absorber extends the domain of control parameters 
that manage a stable periodic radiation processes. An external force suppresses self-organization 
processes. A combined influence of both external force and nonlinear absorber results to more 
complicated picture of self-organization with two reentrant Hopf bifurcations. 
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I. INTRODUCTION 

The most intriguing phenomena in systems 
with nonlinear dynamics is a transition to the 
regime with dissipative structures formation. A 
related problem of such effects investigation in 
systems with large numbers of freedom degrees 
attracts an increasing attention in last three 
decades. Due to self-organization effects a num- 
ber of freedom degrees is reduced and descrip- 
tion of the system dynamics can be performed 
in terms of macroscopic variables. A typical 
picture is realized in laser systems where de- 
scription is provided with a help of amplitudes 
of electric field (or intensity of the radiation), 
polarization and population inversion 

Laser theory shows that corresponding dissi- 
pative structures define formation of pulse or 
modulated signals in homogeneous systems or 
spirals in spatially extended ones 0] . In practice 
a formation of stable periodic radiation can be 
induced by introducing an additional medium 
with nonlinear properties which are realized as 
an absorber or modulator. Such type of lasing 
is known as passive one. Usually, such a kind 
of medium leads to nonlinear dependence of the 
relaxation time of the electric field amplitude 
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[3, 0| or nonlinear dependence of refractive ex- 
ponent composite material can be used 
to introduce different type of such nonlinearity 
■ A coherent dynamics of two- level laser 
systems in the presence of dispersive and ab- 
sorptive effects was observed theoretically and 
experimentally [H, Statistical properties 
of self-organization effects of such type systems 
was discussed in JJl] . Regimes of optical para- 
metric oscillation in a semiconductor microcav- 
ity are studied in . It was found that station- 
ary behaviour of polarization can be described 
by the formalism of non-equilibrium transitions, 
where bistability is observed (see It was 
shown that an oscillating lasing is realized inside 
a bounded domain of the system parameters. 
Another (active) way to initiate a coherent las- 
ing is an introducing an external influence on 
nonlinear processes in the cavity An ac- 
tual problem in laser physics is to find possible 
mechanisms and to set a range of control param- 
eters that manage properties of stable periodic 
radiation (see for instance [1, |3| and citations 
therein). Despite this problem is still opened 
in deterministic (regular) systems a lot of at- 
tention is paid to find coherent regimes under 
influence of stochastic sources [2l| - |24| . 

In this Paper, we are aimed to investigate the 
dynamics of the solid-state class-B laser sys- 
tems which are simple in realization and are 
wide used in physical applications. We consider 
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deterministic models only. According to the 
theoretical approach, based on Floquet analy- 
sis, we will explore in what a manner a non- 
linear medium can induce a stable periodic ra- 
diation. It will be shown that varying in a 
saturation amplitude of electric field and ab- 
sorption coefficient one can arrive at stable and 
unstable dissipative structures. Properties of 
self-organization process induced by an exter- 
nal force influence will be considered. At last, 
a combined influence of both nonlinear medium 
and external force on the system dynamics will 
be described. 

The paper is organized in the following man- 
ner. In Section II we present a model of our 
system where we introduce theoretical construc- 
tions to model an influence of both an absorber 
and external force. Section III is devoted to de- 
velopment of the analytical approach to study 
process of dissipative structures formation. In 
Section IV we apply the derived formalism to in- 
vestigate properties of stable periodic radiation 
in the presence of the absorber, external force 
and its combined affect. Main results and per- 
spectives are collected in the Conclusion (Sec- 
tion V). 

II. MODEL 

Considering a prototype model for a two-level 
laser system, one deals with dimensionless vari- 
ables such as: an electric field amplitude E, 
polarization P and 5 to be a population in- 
version. A standard technique usage Q al- 
lows to reproduce evolution equations for these 
three macroscopic freedom degrees from both 
the Maxwell-type equation for electro-magnetic 
field and density matrix evolution equation. It 
leads to the system of Maxwell-Bloch type that 
is reduced to the Lorenz-Haken model in the 
form 

>c-^E = -E + P, 
ll^P = -P + ES, (1) 
7|^i5 = {S, -S)- EP 

For the single mode laser system a relaxation of 
electric field amplitude E is addressed to losses 



in a bulk of the medium and characterized by 
the velocity >c = 1/2tc, where Tc is a life-time of 
a photon in a cavity. 7j_ is a relaxation velocity 
of nondiagonal elements of density matrix which 
is related to the half-width of a spectral line. 
The relaxation scale for the population inver- 
sion is determined by the velocity 7|| defined by 
both transition probability between two energy 
levels and a corresponding frequency. Se con- 
trols the pump intensity, as usual. The model 
p]) shows a linear combination of the amplitude 
E and polarization P, despite the evolution of 
both P and the pump intensity S are nonlin- 
ear. It is principally that the positive feedback 
of E and S leads to instability in the polariza- 
tion that induces a self-organization. According 
to the Le-Shatelier principle such positive feed- 
back is compensated through negative one in 
third equation (the last term). 

To make an analysis we pass to dimensionless 
variables r = t>c, a = >?/7± and e = 
Hence, the system ([IJ takes the form 

E = -E + P, 

aP = -P + ES, (2) 
eS = {Se - 5) - EP. 

Assuming different combinations between relax- 
ation scales a and e, one can describe three 
possible classes of laser systems. At e,a ^ 1 
we arrive at the laser models of class-A (or- 
ganic dye lasers) with one-dimensional phase 
space, where systems states are represented 
by fixed points only. Here self-organization 
effects are described by a formalism of non- 
equilibrium phase transitions. Class-B (solid- 
state lasers) is characterized by a condition 
CT <C 1. Here phase space is two-dimensional 
and transition processes are of oscillation type 
and hence self-organization processes result in 
dissipative structures formation. For the class- 
C (molecular gas lasers) we set cr, e ^ 1 and 
in three-dimensional phase space a strange at- 
tractor can be realized. At last, the class- 
D (beam masers) is characterized by condition 
cr, e » 1. In this Paper we consider the class-B 
only, where the polarization P is assumed to be 
a microscopic quantity and should be treated 
as fast variable which follows the electric field 
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amplitude E evolution. Such situation is re- 
alized in single mode solid laser systems with 
low-doped crystals (^^203 : Cr^~^) and glasses 
(soda-lime glass), some gas lasers (CO2), fiber 
and semi-conductor lasers [1, 0, [2^ . 

Assuming conditions j± ^ >c, 7|| , one can 
use the adiabatic elimination procedure which 
yields the relation P — ES. As a result, instead 
of the system ([2]) we obtain a two-component 
model in the form 



E = -E{l-S), 

S = £-1 [Se - 5(1 + E^)] 



(3) 



The model ([3|) can not show the stable os- 
cillating regime of the electric field E, itself. 
It was shown experimentally and theoretically 
[li 0j 0| that stable oscillations can be realized 
if an additional nonlinear medium is introduced 
into the cavity. The first way to get the periodic 
lasing is to use a passive modulating medium 
(nonlinear material) to absorb a weak radiation 
and transmit signal with large amplitude. Such 
a type of absorbers is realized in practice as ph- 
thalocyanine fluid in Fabry-Perot cavities |26| . 
gases SFe, BaCh and CO2 [13, Hi. To de- 
scribe action of the absorber it was proposed to 
introduce a nonlinear damping into evolution 
equation for the electric field [29j 



Ik 



kE 



1 + E^/E^' 



(4) 



here Eg is the saturation amplitude. The second 
way is to use an additional medium with nonlin- 
ear refractive exponent n = n [E) 1,01. Such 
type of modulator can be used to increase the 
Q-f actor of laser. We will model action of such 
an effective medium by the external force /e {E) 
assumed in the form 



f.^-A-CE'', 



(5) 



that correspond to action of a bare potential 
V = AE + CE^/3, where coefhcients A, C con- 
trols photon processes in the modulator. We 
use an general construction ^ in order to in- 
vestigate an influence of parameters A and C 
on lasing. In physical applications one can as- 
sociate —A as incident field amplitude, C can 



control nonlinear properties of the refractive in- 
dex n{E). One of the simplest situations is con- 
sidered in [1^, where only a case of A < was 
investigated. 

Combining all above suppositions into the one 
model for a single-mode laser system, we will get 
the generalized system of nonlinear equations 
type of 



E : 

5* = 



UE), 



(6) 



= -Eil -S) + fe{E) 4 
£-1 [Se-Sa + E^)] 

Using two type of additional medium in the cav- 
ity, one can expect that some combinations of 
parameters for both modulator and absorber 
should exist to provide the stable periodic ra- 
diation of the laser. 



TIL MAIN EQUATIONS 

To find mechanisms which takes care of the 
stable dissipative structures formation we will 
use the standard procedure to analyze condi- 
tions where bifurcation into limit cycle occurs 
[30| . To this end we rewrite the system © in a 
most general form 



= /(!)(£;, 5), 

eS^f(^\E,S), 
where effective forces are as follows: 



(7) 



kE 



/(I) {E, S)^-[A + CE^] -E--^ ^,^^2 

/(2)(^,5) =£-1 [Se-Sil + E^) 



ES. 



(8) 

here constructions (H]), ^ are used. 

Wc deal with a problem of nonlinear dynam- 
ics and present a behaviour of the system in the 
phase plane {E, S) . Firstly, we consider steady 
states Eq and Sq , defined as coordinates of fixed 
points in the phase plane. Setting _E — and 
S = 0, one can find steady states as solutions of 
stationary equations 

^0 ( — ^ - J^L. ~ '2CEo -l]^A. 



1 



El 



El 



E^ 



So — <5'e(l + Eq 



(9) 
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A behaviour of phase trajectories in the vicinity 
of these fixed points can be analyzed with a help 
of the Lyapunov exponents approach. Here time 
dependent solutions of above system are as- 
sumed to be in the form E oc e^*, A = A + ioj, 
where A controls the stability of the phase tra- 
jectories, u; determines pulse frequency of the 
signal. Magnitudes for real and imaginary parts 
of A are calculated according to the Jacobi ma- 
trix elements 

M,, = (^) ; = {E, S}, t,j = 1, 

(10) 

where subscript relates to steady states. In- 
serting ([8]) into definition (fTO]) . we get matrix 
elements 



Mil = -Mo + ^o, 
Mn = 1 + 2CEa + K- 



(11) 



\l + El/E^,Y' 
Mi2 = So; M21 = -2e-^SoEo; (12) 
M22^-e-\l + E^). 

Then, an equation for eigenvalues and eigenvec- 
tors 



gives expressions for A and loq as follows: 

A=^[(So-Afo)-e-Ml + ^o)], 



(13) 



To investigate a stability of such a limit cy- 
cle we analyze a behaviour of trajectories in the 
vicinity of the fixed point {Eg, Sq). To this end 
we rewrite motion equations ([7]) where variables 
E and S are count off from stationary magni- 
tudes Eo,So- To do this one can use following 
transformation 

X = Xo + P-S, (17) 
where notations for pseudovectors are used: 

S = 



c^o = -]j8e-^S„E^ - [(5o - Mo) + e-^l + El)f 

(14) 



If the real part of the Lyapunov exponent A = 
then a fixed point (Eq^Sq) is addressed to a 
center of a limit cycle. It leads to relation 



e{So-Mo) >1 + El, 



(15) 



and yields a condition for the frequency of os- 
cillations 

SeSoE^ > [e(5o - Mo) + (1 + E^)] ' . (16) 



X = 



E — Eq 

S — Sq 



(18) 



The corresponding transformation matrix P is 
obtained with a help of eigenvector V compo- 
nents, i.e.: 



P 



V 



Vi 
V2 



(19) 



Assuming Vi = 1, for the second component 
V2 from Eq. p^ one gets 



V2 



(Mo - -So) + itoc 



Eo 



LOc = ^^oIa=o 



2SeEQ 



a+E^or 



-,1/2 



(20) 



Hence, the transformation matrix (jT9|) takes the 
form 



P = 



1 







{Mo - So)/Eo - ujc/Eo 



(21) 



It leads to evolution equations for deviations 
written in a vector form 

= F = p-^f. (22) 

Here a pseudovector of the canonical force 



F 



^(2) 



satisfies conditions |30|-|3 



dF 
85 







(23) 



(24) 
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and has following components: 



(2). 



Ma ~ Sq 



(25) 



(26) 



Above procedure allows to find the stability of 
the manifold formed by the fixed point (i?o, Sq). 
Using the standard technique [3(]|, one can say 
that the limit cycle is stable only if a real part 
of the Floquet exponent 



gl502| 



1 

2521, 



* ^ 2^ (^ffii320 - 2|5ii| 

(27) 

is negative in a bifurcation point. Structure 
constants in the definition ((27|) are described 
by derivatives with respect to E and S*, denoted 
with subscripts: 



511 ^ 1 



502 
520 



(1) 



EE 



5^(1) 
SS 



(2) 



^ EE ^ ' 



EE 



SS + ^^BS 



5.(2) 



'( 



EE 



^ SS ^ '^^ ES 



(28) 



(29) 



p(2) 
SSS 



521 - g I {FeIe + ^BSs) + i^^EES + 

i {Peee + Pess) ~ (^fiis + P'sss) } ■ 

(30) 

Using some algebra, the stability condition for 
the limit cycle can be written as follows 

2a{iP^ - Cf + al3eSa{l + 2/3£So) + + I3e) 

[C ~ V'«)(a2 - 1 + 2(3eSQ + 2af3eEo), 

(31) 

where notations 

kEs^ E {~3 E,^ + E^) 



IV. ANALYSIS OF HOPF BIFURCATIONS 
A. Influence of nonlinear absorber 

To proceed let us consider steady states be- 
haviour under supposition that action of the ab- 
sorber is given by expression (|4|), /e = 0. Set- 
ting £' = 5 = 0, one gets stationary values of 
the electric field amplitude Eq shown in Fig[TJ 
A steady states analysis allows to find that a 
bistable regime is realized only if k < Kmin , here 
Kmm = E1/{1 — E1). In such a case one gets the 
hysteresis loop in Eo{Se) dependence in the do- 
main [S'ccS'c] (curve 1) which disappears when 
the threshold Kmin is crossed, where 



S'c = 1 + K, S'cO = 1 + Es 



(32) 



1-Er 

The behaviour of the amplitude Eq is the same 




< 



6- 



{Es' + E^y 

kEs^ (-6 E^E.^ + E'^ + 



E, 



E^y 



are used. 



FIG. 1: Stationary amplitude -Eo vs. pump inten- 
sity Se at different values of absorption coefficient 
at Es = 0.9: curve 1 - k = 10.0; curve 2 - k = 4.0 

as in the first order phase transitions where zero 
value of Eq below S'cO corresponds to a disor- 
dered state, values £0 7^ (solid line) relate 
to an ordered state, whereas intermediate mag- 
nitudes of Eo (dotted line) correspond to un- 
stable state. The critical value for the absorp- 
tion coefficient is realized only if the saturation 
amplitude Eg < 1. In opposite case one can 
get the stationary picture of the second order 
phase transition where Eq increases monoton- 
ically from if the critical value Sc is crossed 
(curve 2). 
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The analysis of the Floquet exponent aUows 
to find the phase diagram (Fig[2l), which shows 
the stable periodic radiation (formation of limit 
cycle in the phase plane {E, S)). In FigElthe do- 




FIG. 2: Phase diagram in the presence of the non- 
linear absorber at Es =0.9 



main I defines configuration of the phase space 
with both a stable focus (ordered state) and a 
saddle point (disordered state); in the domain II 
only disordered state is realized (node); the do- 
main III is characterized by the hysteresis loop, 
where ordered state corresponds to unstable fo- 
cus, unstable state is represented by a saddle, 
disordered state is a node. Inside the domain IV 
the stable limit cycle is formed (Fig[3^), which 
transforms into stable focus, unstable and sta- 
ble cycles if dotted line is crossed (FigOjD). An 
influence of the parameters of the absorber on a 
topology of phase plane is shown in FiglH Here 
an increase in the absorption coefficient k at 
small Eg leads to transformation of unstable fo- 
cus into a stable one with additional node and 
saddle points appearing. At values k and Es, 
corresponding to the dashed line, one gets an 
unstable limit cycle and in the domain bounded 
by dashed and solid lines one gets the unstable 
focus, node and saddle. When the solid line is 
crossed the phase portrait is characterized by 





b) " 



FIG. 3: Dissipative structures corresponding to do- 
main IV (a) and dotted line (b) in Fig[2] Stable 
and unstable limit cycles: (a) — Es = 0.9, k = 8.0, 
Se = 9.9; (b) Es = 0.9, k = 10.0, Se = 12.3 
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FIG. 4: Transformation of topology phase plane un- 
der influence of the absorber parameters k. Eg at 
Se = 5.0 



a single node. An increase in k at saturation 
amplitude -Es ~ 1 transforms an unstable fo- 
cus into a stable limit cycle, which becomes un- 
stable at values that correspond to the dashed 
line. In the domain bounded by the dashed and 
straight horizontal lines there is a single unsta- 
ble focus only. A further increase in k trans- 
forms this focus into a node. 

The frequency of pulse radiation regime ap- 
pears at non zero value, that correspond to the 
first bifurcation point Se a further increase in 
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the pump intensity, leads to the growth of lOc 
till the second critical point Se is achieved. We 
have analyzed behavior of pulse radiation fre- 
quency at different values of the absorption co- 
efficient K. According to FigIS] an increase in k, 
at fixed saturation amplitude magnitudes leads 
to the shift of minimal and maximal values of 
LOc despite a topology of the dependence LodSe) 
is not changed. Obtained results are in good 
corresponding with exp erimental observations 
of such dependence [34| . 




e 



FIG. 5: Dependence of the stable pulse radiation 
frequency ujc vs. pump intensity Se for the system 
with absorption effect at Es = 0.9 and k = 5, 6, 7, 
8, 9 (from the left to the right). Insertion shows a 
typical stable pulse signal 

Therefore, the dispersion in the relaxation 
time of the electric field amplitude promot- 
ing by the absorbing influence, leads to forma- 
tion of the stable periodic radiation at satura- 
tion amplitude i?s — 1. 
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FIG. 6; Stationary amplitude Eo vs. Ss at A — 1.9, 
C = -0.455 



can be formed only if a stable focus is trans- 
formed into an unstable one and vice versa (see 
Figil). Here at Se < Sc and Se > S" the 
phase portrait is characterized by single sad- 
dle point 5*1 or S2, respectively. In the do- 
main Sc < Se < S^ one gets two saddles Si 
and 5*2, divided by an unstable focus F^- If 
S" < Se < S"^, then such saddles are divided 
by a stable focus Fg. Only if Se — 5*° we will 
get a trivial situation, where K$ = 0. It means 
a formation of nested loops of neutral stabil- 
ity (Fig[7]). Therefore, external force suppresses 
processes of dissipative structure formation. 



C. Combined effect of external modulator and 
nonlinear absorber 



B. Influence of external modulator 

Let US consider an influence of the external 
source fe at = 0. It is principally important 
that the periodic radiation is possible only if pa- 
rameter that controls nonlinear effects C < 0. 
Here stationary behavior of the field Eq versus 
pump intensity Se is shown in Fig[6l Analysis 
of the Floquet exponent shows that limit cycles 



Now we consider an influence of both external 
modulator and nonlinear absorber on the pro- 
cesses of dissipative structure formation. Set- 
ting E — S — in the system (O, one gets 
stationary values of the electric field amplitude 
Eq shown in Fig[5] As it is seen, if the modu- 
lator is turned oS {A ~ C = 0) then we have 
a single stable state with no radiation at small 
values of the pump parameter Se- If the thresh- 
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FIG. 7: Typical phase portrait of the system with 
external source at Se = 5.19152 A = 1.9, C = 
-0.33025 
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FIG. 8: Stationary value of electric field amplitude 
Eq vs. pumping intensity Se for the system ((Tjl at 
C = 0, K = 5.0 and different magnitudes for the 
parameter of spontaneous emission A 



old given by expression = 1 + k is crossed, 



then a new solution of the steady state equa- 
tion appears and we have a stationary radia- 
tion with an amplitude Eq ^ which increases 
with an increase in the pump intensity. If we 
set A < at C = then we will get a single 
stable solution on the whole axis of the pump 
parameter magnitudes which defines the radi- 
ation amplitude Eq. In the opposite case of 
^ > one gets two stationary solutions, only 
if the energy barrier S^, given by the solution 
of equation = /(A, Kg, k), is overcame. 

Next, we investigate conditions where stable 
periodic radiation can be realized. To this end 
we need to determine a domain defined by con- 
ditions A = and 5R<i> < where periodic solu- 
tions of the system ([7]) are exist. Correspond- 
ing solutions of the Eq. (PT|) are shown in FiglHl 
It illustrates domains of the absorption coeffi- 
cient K and pump intensity Se magnitudes at 
different intensities A, C where the stable radi- 
ation process is realized. As FiglS] shows, if we 
set an absorber inside the cavity only, then a 
semi-limited domain of k and Se magnitudes is 
formed; inside of this domain the stable periodic 
radiation is possible. Introducing a modulator 
with A > 0, C — (see Fig|9^), such a domain 
becomes totally limited. Moreover, an increase 
in the parameter A leads to restriction of the 
values for the collective parameter k and pump 
intensity 5*6, at which one has stable periodic 
radiation. At large values A such domain is de- 
generated into the line. From Fig[5]D one can 
see that an increase in the C at A ~ leads to 
extension of the domain of stable periodic radi- 
ation that occurs at large magnitudes of pump 
intensity parameter. 

An influence of nonlinear processes in the 
modulator on a picture of the stable periodic 
radiation formation is presented in Fig llOl It is 
seen, ii A < then there is only stable station- 
ary state (see Fig|S]) which is a focus (5RA < 0, 
$JA ^ 0) on a phase plane {E, S). Such a fixed 
point is transformed into a manifold if control 
parameters are in the domain including its bor- 
der shown in FigfTUk. Such a manifold is a limit 
cycle (3?$ < 0, = 0) in the phase plane 
{E,S), that attracts all phase trajectories in 
the vicinity of it. From a physical viewpoint 



9 








-0.20 


sf 


0.15 


7.5- 


A 




0.10 




0.05 





FIG. 9: Phase diagrams of periodic radiation ef- 
fect: (a) an influence of the spontaneous emission 
intensity A at C = 0; (b) influence of the photon 
scattering intensity C a,t A — 



it means the formation of the stable pulse pe- 
riodic radiation. The domain shown in Fig llOb 
is limited by the value of intensity of nonlinear 
processes C*. One needs to note that if pho- 
ton scattering occurs with intensities C < 
then an increase in pump intensity Se induces 
formation of stable periodic radiation at magni- 
tude S^^ and destroys it at S^^. In other words, 



a) 




b) 



0,5 



I.O 



FIG. 10: Phase diagrams of Hopf bifurcation: in- 
fluence of the nonlinear processes intensity C at 
K = 5.0: (a) A = -0.1; (b) A = 0.1 



one gets the situation where the only one rea- 
son serves as stimulus for both self-organization 
and desorganization. 

A picture became more complicated at ^ > 0. 
At first let us discuss the phase diagram shown 
in FigllOb. At pump limited by the dashed 
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curve in Fig llOb there are no stationary solu- 
tions and, hence, no stable regimes of radiation. 
Next, processes of spontaneous photon annihi- 
lation reduce a domain of stable periodic radi- 
ation at pump intensities above dashed curve, 
here a domain of unstable behaviour of phase 
trajectories appears. At small C < Co such a 
stationary regime is defined by the correspond- 
ing stationary solution which is an unstable fo- 
cus (5RA > 0, 3 A 7^ 0). The related fixed point 
is defined as an upper branch of the dashed 
curve in Fig [51 At large values C" < C < 
one has a picture similar to discussed above. At 
intermediate values Co < Cr < C" one can get 
a very complicated picture of self-organization. 
Here with an increase in Se we have follow- 
ing picture of transformations: (i) the system 
passes from the unstationary regime to station- 
ary one with fixed point to be a stable focus 
(FigfTTk): (ii) at values Se^^^^ one has a sta- 
ble periodic radiation that exists till magnitudes 
< Sf-^^ (Figdlb); (iii) a further increase m 
pump intensity destroys the limit cycle and a 
system pass to unstable regime which is char- 
acterized by the unstable focus (FigfTTh): (iv) 

if critical value Se^^"^^ is achieved, then a new 
Hopf bifurcation occurs and the system evolves 
according to periodic trajectories fFigfTTb): (v) 
at last, such a coherent regime is destroyed at 

Se > Sf^'^ (FigEJi). 

Let us consider more closely properties of 
phase diagram (FigH]), which shows magnitudes 
of the absorption coefficient k and the control 
parameter Se- Here the thin solid curve (bi- 
furcation line) defines critical magnitudes for k 
and 5*6 where stationary states appeared. In 
the domain with 5RA > one has an unstable 
focus (see FigfTSk'). The stable limit cycle is re- 
alized inside the bounded domain with 5R$ < 
(FigfT3b). At small k and large Se one has a sta- 
ble focus fFigfTHb). Therefore, one gets a trans- 
formation of topology of attractors in the phase 
plane if parameters Se or k are changed. An in- 
crease in absorption coefficient k at fixed pump 
intensity Se will produce the oscillating regime 
(transition from a stable focus to limit cycle). 
Such stable periodic regime can be destroyed at 



large magnitudes of k (transition from the limit 
cycle into repeller — unstable focus) and a fur- 
ther increase in k leads to the absence of any 
stationary regime at all. However, the stable 
periodic solution is observed not on a whole bor- 
der of the indicated domain. Figure [T^] shows 
that a stable dissipative structure is formed in- 
side the domain and on the thick solid lines only 
(5R$ < 0). A part of the domain border plotted 
as dashed line corresponds to conditions A = 
and 9fi<I> > 0, which mean existence of unstable 
periodic solution (see FigfT5h). Hence, there is 
a point where 51?$ — and periodic solution 
changes its stability. In this point the phase 
portrait of the system is characterized by a set 
of nested loops. 



V. CONCLUSIONS 

In this Paper we have analyzed properties 
of self-organization processes in the two-level 
class-B laser systems in the presence of absorp- 
tion effects and influence of the external force. 
We have shown that due to the nonlinear damp- 
ing the domain of control parameters of the cav- 
ity with the stable pulse radiation is realized. It 
was shown that varying a saturation amplitude 
and absorption coefhcient one can pass to dif- 
ferent type of radiation, characterized by fixed 
point in the phase space type of: stable and un- 
stable focuses, stable and unstable limit cycles. 
Introducing the external force that leads to ad- 
ditional nonlinear effects that reduce domains 
of control parameters with stable periodic radi- 
ation. It is principally important that due to 
the external force influence one can get reen- 
trant Hopf bifurcation. Here there is a wide 
range of the external force parameters, where 
both stable and unstable dissipative structures 
are in the phase space. Our results are in good 
correspondence with theoretical ones yJljl and 
experimental observations [1, [13, H^, l34l. |35|| . 

In our investigation we have considered the 
simplest case, where relaxation velocities of the 
electric field and population inversion are of the 
same order. In real systems of the solid-state 
class-B lasers e = ~ 10~^ 10~^, in 



11 



gas lasers of such class e ~ 1. As was shown 
theoretically and experimentally a differ- 
ence between above relaxation velocities will 
not change the picture of stable pulse regime 
qualitatively. Experimental investigation shows 
quantitative changes only. 

In our consideration the construction for the 
external force can be applied to describe influ- 



ence of the nonlinear processes: in the nonlinear 
medium with the nonlinear dependence of re- 
fractive index (a variation of the parameter C) ; 
introducing an external incident field with am- 
plitude A < 0; more complicated picture with 
arbitrary A and C under supposition of the dy- 
namic system stability only. 
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FIG. 11: Phase portraits related to the domains in 
FigdOjj: (a,)- A = 0.1, C = 0.5, k = 5.0, Se = 8.8; 
(h) - A = 0.1, C = 0.5, K = 5.0, Se = 7.5; (c) - 
A = 0.1, C = 0.5, K = 5.0, 5"^ = 7.0 
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FIG. 12: Phase diagram at > 0, C = 



d) 



FIG. 13: Phase portraits related to the phase dia- 
gram in Figdal (a) - A = 0.1, C = 0.0, k = 5.1, 
Se = 8.0; (b) - A = 0.1, C = 0.0, k = 5.5, Se = 8.0; 
(c) - A = 0.1, C = 0.0, «: = 5.6, & = 8.0; (d) - 
A = 0.1, C = 0.0, K = 9.86, = 14.0. 



